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SECTION  I 


INTRODUCTION 

The  most  frequently  encountered  crack  geometry  in  the  aerospace 
field  is  the  surface  flaw  in  bodies  of  finite  thickness.  These  flaws 
may  be  initiated  in  a variety  of  ways  from  surface  imperfections,  and 
their  growth  under  cyclic  loading  is  believed  to  result  in  more  failures 
of  aerospace  structural  components  than  from  any  other  cause.  Although 
an  analytical  treatment  for  the  local  field  around  a surface  flaw  in  a 
half  space  was  provided  by  Irwin  in  1962  (1)  the  problem  of  the  surface 
flaw  in  a body  of  finite  thickness  has  resisted  efforts  of  applied 
mathematicians  to  render  it  tractable  in  the  full  field  closed  form 
sense,  even  for  the  elastic  case.  In  recent  years,  capitalizing  on  the 
speed  and  storage  capacity  of  modern  digital  computers,  various  numerical 
techniques  (2)^(3)  have  been  focused  on  the  problem  and  progress  in  this 
direction  is  continuing  (4).  However,  rigorous  convergence  proofs  for 
three  dimensional  cracked  body  problems  are  difficult  to  achieve,  and 
research  programs  are  now  uniformly  requiring  computer  code  verification 
on  an  experimental  basis  in  order  to  validate  numerical  analyses. 

Several  years  ago,  the  senior  author  and  his  associates  began  to 
develop  (5-7)  an  experimental  technique  based  upon  an  idea  of  G.  R. 

Irwin  for  direct  measurement  of  stress  intensity  factors  (SIFs)  photo- 
elastically.  They  adapted  the  idea  to  the  so  called  "frozen  stress" 
method,  which  they  then  combined  with  an  elementary  least  squares 
digital  computer  program  for  estimating  SIF  distributions  for  three 
dimensional  cracked  body  problems.  The  early  form  of  the  technique  was 
identified  by  the  Society  for  Experimental  Stress  Analysis  as  a viable 
approach  to  three  dimensional  cracked  body  problems  and  was  featured  in 
the  lead  chapter  of  a recent  SESA  monograph  (8)  on  Experimental  Fracture 
Mechanics.  The  method  has  been  modified  and  refined  since  that  time  (9- 
15)  and  has  been  successfully  applied  to  several  surface  flaw  problems 
to  date  (7),  (16-30). 

A problem  which  has  received  only  limited  attention  (31)  in  the 
literature  deals  with  the  surface  flaw  in  a finite  body  under  mixed  mode 
loading.  Absence  of  information  on  this  subject  is  not  unexpected  since 
this  problem  is  substantially  more  complex  than  the  pure  Mode  I problem 
discussed  above.  The  study  described  in  the  sequel  is  concerned  with 
the  application  of  the  frozen  stress  technique  noted  above  to  the  Mixed 
Mode  Surface  Flaw  problem  in  which  part-circular  flaw  geometries  were 
employed. 
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SECTION  II 


ANALYTICAL  CONSIDERATIONS 


The  main  quantity  obtained  from  phutoelastic  analysis  is  the 
maximum  shearing  stress  in  a plane  normal  to  the  direction  of  travel  of 
the  polarized  light  field.  In  this  section  we  develop  the  analytical 
background  for  converting  this  information  into  SIF  values. 

It  has  been  shown  (32)  that  the  stresses  near  the  border  of  an 
elliptical  crack,  when  expressed  in  coordinates  in  a plane  perpendicular 
to  the  border,  have  the  same  form  as  the  stresses  in  a plane  perpen- 
dicular to  the  border  of  a crack  with  a straight  front.  For  a crack 
with  a circular  front,  the  same  type  of  results  had  been  obtained  by 
Sneddon  (33).  This  form  of  the  stresses  is  further  verified  by  the  work 
of  Thresher  and  Smith  (34)  for  the  part  circular  crack. 


Consider  an  infinite  plate  with  finite  thickness  which  is  subjected 
to  uniaxial  loading.  The  stress  distribution  near  a part-circular 
crack,  which  is  inclined  to  the  plate  boundary  at  an  angle  3,  can  be 
found  using  the  following  arguments. 

The  geometry  is  shown  in  Figure  1.  Based  on  previous  experimental 
work  the  infinite  plate  assumption  appears  valid  for  2m  < 0.30  w and 
L > 3w  where  2m  is  the  crack  length,  w is  the  plate  width,  and  L is  the 
plate  length. 

The  orthogonal  coordinate  system  nzt  used  in  this  problem  (see 
Figure  2)  is  a local  frame  and  it  will  move  along  the  leading  edge  of 
the  part  circular  crack,  such  that  z will  always  be  perpendicular  to  the 
crack  surface,  and  n and  t will  define  the  crack  surface  and  will 
remain  perpendicular  and  tangent  respectively  to  the  flaw  border. 

It  is  clear  that  the  three  modes  of  loading  Mode  I,  Mode  II,  and 
Mode  III  are  present  at  any  point  along  the  crack  front. 


The  stress  distribution  near  the  part-circular  crack, 
to  the  opening  mode  of  deformation  can  be  taken  as  (15) . 
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where  K1  is  the  Mode  I SIF  and  the  coordinates  r and  0 are  shown  in 
Figure  2.  ol  represent  the  contribution  of  the  Mode  I regular  stress 
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field  to  the  measurement  zone  which  is  taken  far  enough  from  the  crack 
tip  to  avoid  a non-linear  zone  very  near  the  tip.  While  they  may 
generally  be  regarded  as  expressible  in  Taylor  Series  Expansions,  it 
turns  out  that  only  the  leading  terms  of  said  series  are  necessary  so 
that  are  constants  for  a given  point  in  the  flaw  border  but  vary 
from  point  to  point.  Stresses  corresponding  to  the  Mode  II  can  be  t 

K 


as 


nn 


zz 


nz 


taken 


(27Tr) 


. 0 ro  . 9 30  2 

YJ2  sln  y L ^ + cos  j cos  — jj-a^ 


(27rr) 

K0 


. 0 0 30  2 

rrj ~ sm  tt  cos  tt  cos  — -a 
1/2  2 2 2 zz 


0 r i , 0 . 30n  2 

T-J7T  COS  TT  l 1 - Sin  7T  S±U  7? J — T 

1/2  2 2 2 nz 


(2) 


(27Tr) 

2 1 
where  K2  and  j are  analogous  to  and  a^j 


And  finally  the  stresses  corresponding  to  the  Mode  III  loading  can 
be  taken  as 
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with  and  analogous  to  K-^  and  a^j 


If  the  above  modes  of  loading  are  superimposed,  one  gets  the  fol- 
lowing stress  distribution  in  a plane  perpendicular  to  the  crack  front 
(the  n-z  plane) 
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which  are  independent  of  Equations  (3). 

While  has  no  influence  upon  the  singular  stress  field  itself, 

it  does  alter  the  isochromatic  fringe  pattern,  which  is  proportional  to 
the  maximum  in-plane  shearing  stress. 


From  the  stress  field  given  in  Equations  4,  the  maximum  shearing 
stress  in  the  plane  perpendicular  to  the  crack  front,  nz , can  be  ob- 
tained using 
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and,  truncating  to  the  same  order  as  Equation  4,  one  gets 
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and 


B = B (0,  Kr  K2,  0°  ) 


Photoelastically , the  maximum  shearing  stress  (in  the  nz  plane)  is 
determined  from  the  stress  optic  law 


fnT 

T = TZT 
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where 


f is  the  photoelastic  fringe  constant  for  the  material 
t*  is  the  thickness  of  the  slice  measured  parallel  to  the 
direction  of  the  light,  and 
nf  is  the  isochromatic  fringe  order. 

Now,  in  general,  the  effect  of  involves  both  a folding  and  a 

change  in  eccentricity  of  the  fringe  loops  (13).  If  folding  occurs,  0m, 
the  angle  along  which  the  distance  to  a fringe  from  the  crack  tip  is 
greatest,  will  vary  with  the  fringe  order  nT  and  one  must  plot  0m  vs 
r/a  and  extrapolate  to  the  origin  in  order  to  obtain  0°,  the  value  of  0m 
associated  with  and  K^.  In  the  present  problem  0m  was  constant  over 
the  data  range  in  the  fashion  indicated  qualitatively  by  Figure  3. 
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Upon  computing 
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Since  0m  can  be  measured  experimentally  (I^/K-^)  can  be  calculated  from 
Equation  9.  Then  by  combining  the  stress-optic  law  with  a modified  form 
of  Equation  6. 
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and  one  can  solve  for  the  individual  values  of  and  K^. 


In  order  to  do  this,  one  must  obtain 

K*  = [ (K,  sin  0°  + 2K0  cos  0°)2  + (K0  sin  0°)2]1/2 
1 m z m z m 


from  by  plotting  K*^  = (Tnz)ma  (Sttt ) vs  (“)  identifying  a 
linear  zone,  and  extrapolatingZtoa£he  origin,  which  will  yield  K*.  A 
typical  set  of  fringe  data  and  SIF  determination  is  given  in  Figure  4. 
Once  K*,  K^/K^,  and  0^  are  known  and  can  be  calculated  which  then 
can  be  normalized  using  proper  quantities. 


Note  that  the  above  approach  utilizes  a two  parameter  (A,  B) 
model,  since  the  linear  zone  of  Equation  10  can  be  located  experimen- 
tally. However,  if  one  cannot  locate  such  a zone  experimentally  then 
additional  terms  leading  to  an  equation  of  the  form 

oo 

Tmax  1/2  + ^ B r11^  (12) 

r n=o  n 


with  suitable  truncation  criteria  must  be  considered.  Since  such 
criteria  are  not  yet  established,  this  latter  approach  is  avoided  where 
possible  and  was  not  necessary  in  the  studies  described  in  the  sequel. 

The  stress  distribution  in  a plane  perpendicular  to  the  crack  sur- 
face and  tangent  to  the  crack  front  (zt  plane)  or  in  a plane  parallel  to 
the  zt  plane  can  be  found  from  Equations  4. 
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In  order  to  arrive  at  a value  for  azz,  prior  experiments  by  the 
authors  indicate  that  the  usual  assumption  of  plane  strain  (for  the 
plane  problem)  may  not  be  valid  here.  However,  if  one  assumes  a state 
of  nearly  generalized  plane  strain  such  that  the  value  of  ett  can  be 
considered  constant  over  a portion  of  the  length  of  the  flaw  border, 
then  the  observed  state  of  varying  transverse  constraint  along  the  flaw 
border  can  be  approximated  rather  well.  Thus,  we  assume 
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flaw  border  and  where 


a 

nn 


K., 


(27rr) 


0 ri  . 0 . 36, 

YJ2  cos  ~2  L1  “ sin  *2  Sln  ~2" ' 


2 . 0 T9  a.  0 301 

sm  y [2  + cos  j cos  — J - a 


(27Tr) 


1/2 


o 

nn 


For  v = y*  (as  in  the  present  experiments) 
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Moreover,  from  Mode  III,  we  have, 
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Consider  the  line  normal  to  the  crack  surf ace^which  passes  through 
the  crack  tip  in  the  zt  plane.  For  this  case  0 = y,  and  when  substi- 
tuted in  Equations  13,  13a  and  13b,  there  results 
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then 


^Tzt^max  1/2 
r 


+ D 


(15) 


where 


c - —TIT!  'T5  (Ki  + k2)2  + k3J1/2 

(4tt) 


and 


D = D (K1,  K2,  K3,  Ee  and  o°  ) 

Now  by  combining  the  stress  optic  law  with  a modified  form  of 
Equation  15 


(t  ) 

zt  max 


fn' 
2t 1 


K** 

AP 


(87Tr) 


1/2 


where 


KAP  ■ ^ [li  (K1AP  + K2AP)2  + K3AP]1/2 


the  value  of  K3  can  be  obtained. 

In  order  to  do  this,  one  must  obtain 

K**  = /2  {y|-  (Kj.  + K2)2  + K2}1//2  from  K**  by  plotting  a 

K**  = (t  ) (Sur)1^2  vs  (r/a)1^2  , 

AP  zt  max 

identifying  a linear  zone,  and  extrapolating  to  the  origin,  which  will 
yield  K**. 


7 


SECTION  III 


THE  EXPERIMENTS 

"Frozen  stress"  experiments  require  that  the  model  material  exhibit 
certain  special  characteristics  in  addition  to  the  optical  transparency 
required  of  all  photoelastic  materials.  Stress  freezing  materials 
behave  in  a Kelvin-like  manner  (Figure  5)  at  room  temperature.  However, 
when  heated  to  a certain  temperature,  known  as  the  "critical"  tempera- 
ture, the  anelastic  coefficient  y-K),  and  the  material  becomes  linearly 
elastic,  exhibiting  an  elastic  modulus  over  two  orders  of  magnitude 
below  the  room  temperature  value  and  an  isochromatic  sensitivity  twenty- 
five  times  greater  than  at  room  temperature.  If  the  cracked  body  is 
heated  to  critical  temperature,  loaded  and  cooled  slowly  under  load,  the 
recovery  effects  upon  unloading  are  negligible  due  to  the  difference  in 
modulus  and  isochromatic  sensitivity  between  critical  and  room  tempera- 
ture. In  fact  the  model  may  be  sliced  (and  this  is  done  parallel  to  the 
nz  plane  [Figure  1])  without  altering  the  isochromatic  or  deformation 
patterns  imposed  above  critical  temperature.  Thus  "frozen  stress"  is  a 
misnomer  in  the  sense  that  live  or  macrostresses  are  not  trapped  in  the 
body,  but  rather  deformation  fields  and  the  isochromatic  evidence  of 
live  stress  above  critical  temperature. 

A series  of  19  photoelastic  models  were  prepared  in  accordance  with 
Table  I.  Cracks  were  simulated  with  thin  sawn  slits  tapered  with  a 30° 
vee  tip  to  a root  radius  of  2.5  x lCT^mm.  Such  geometries  had  been 
shown  analytically  (35)  and  experimentally  (10)  to  yield  SIF  estimates 
to  within  2%  of  natural  cracks.  Mixed  mode  fields  were  then  frozen  into 
17  of  these  models  and  primary  slices  were  removed  near  the  crack  tip 
along  the  flaw  border  parallel  to  the  nz  plane.  Isochromatics  were  read 
from  these  slices  by  placing  them  in  a circularly  polarized  white  light 
field  in  a bath  of  matching  index  fluid,  and  reading  tint  of  passage 
using  the  Tardy  Method  along  predetermined  directions.  Typical  fringe 
patterns  are  shown  in  Figure  6.  Optical  data  were  then  fed  into  a least 
squares  digital  computer  program  for  establishing  the  curve  of  best  fit, 
from  which  K*  estimates  were  obtained  by  extrapolation  as  in  Figure  4. 
These  values  were  then  converted  to  and  K2  values  analytically. 
Secondary  slices  were  then  obtained  from  the  primary  slices  and  analyzed 
in  a similar  manner  for  K3. 
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SECTION  IV 


RESULTS 

Test  geometries  and  loads  for  all  sixteen  experiments  are  given  in 
Table  I according  to  increasing  a/T1.  Test  4 represents  an  average  of 
four  identical  tests  and  so  nineteen  tests  were  conducted  altogether. 

SIF  distributions  for  Modes  I and_II  are  presented  in  Figures  7 through 
22.  The  normalizing  factor  = a (ira) l/^sin^B , is  the  Mode  I SIF  for  a 
through  straight  front  crack  in  an  inf inite_body  whose  plane  is_inclined 
at  an  angle  3 to  the  remote  uniform  stress  a.  Similarly,  K2  = a( Tra)l/2 
sin  3 cos  3 is  the  Mode  II  SIF  for  the  same  problem.  Consequently,  the 
normalizing  factors  do  not  vary  along  the  flaw  border. 

There  are  several  points  which  the  authors  wish  to  make  regarding 
the  results  of  Figures  7 through  22.  First,  it  has  been  pointed  out  by 
other  investigators  (36),  (37)  that  the  calculations  used  in  this  type 
of  approach  are  quite  sensitive  to  values  of  0m,  especially  for  0m 
I 90°  (Fig.  3).  In  order  to  assess  this  effect  upon  the  present  results, 
four  "identical"  tests  were  conducted  (Tests  4)  for  the  geometry  expected 
to  produce  the  highest  scatter  due  to  the  sensitivity  to  0m.  Results 
are  given  in  Figure  7 showing  the  average  values  with  a scatter  band  for 
the  data  from  the  four  tests.  These  results  indicate  that,  indeed,  when 
the  method  is  applied  to  mixed  mode  studies  where  0m  is  measured  rather 
than  assumed  (as  for  Mode  I)  the  resulting  accuracy  in  measuring  both 
Modes  I and  II  is  reduced  from  the  usual  scatter  of  + 5%  for  the  pure 

Mode  I case  to  as  large  as  + 12%  for  Mixed  Mode  problems  with  artificial 

cracks  at  the  point  of  maximum  flaw  penetration. 

Due  to  limitations  on  cutter  sizes  for  producing  the  part  circular 
flaws,  it  was  not  possible  to  hold  two  of  the  variables  3,  a/ Tf,  a/2m 
constant  while  varying  the  other  one.  In  order  to  cover  a range  of 
interest,  three  levels  of  a/Tf  were  treated,  two  values  of  3 and  two 
levels  of  a/2m.  With  this  approach,  it  turned  out  that  when  3 was 
increased  from  30°  to  60°,  an  increase  in  a/Tf  also  occurred.  Although 
the  influence  of  the  inclination  of  the  remote  stress  to  the  crack 
surface  was  normalized  out  using  K values,  the  influence  of  the  inclina- 
tion of  the  plate  surfaces  to  the  Mixed  Mode  Components  remains  embedded 

in  the  data.  In  four  of  the  tests  (13  thru  16)  Figs.  20-22,  a study  of 
the  influence  of  2m/w  was  conducted. 

Within  the  above  limitations,  the  following  trends  can  be  observed: 

Effect  of  a/T T , 3 - Tests  2,  4,  6,  7 and  11  (Figures  9,  7,  12,  13,  18 
respectively)  show  a continuously  increasing  a/T!  with  little  change  in 
a/2m.  Some  of  these  tests  (2  and  6)  were  for  3 = 30°  and  the  others  (4, 

7 and  11)  were  for  3 = 60°.  Viewing  these  tests  collectively,  the  only 
consistent  result  appears  to  be  that  maximum  values  of  K^/K^  and  K2/K2 
do  not  occur  at  maximum  flaw  penetration  for  any  case.  No  significant 
trends,  however,  are  noted.  On  the  other  hand,  if  one  separates  the 
tests  according  to  values  of  3,  (i.e.  2 and  6 and  4,  7,  11)  definite 
trends  do  appear.  In  both  sets,  the  magnitude  of  Ki/K]^  increases  with 
af/T  and  shifts  away  from  the  point  of  maximum  flaw  penetration.  Also, 
while  Ki/K-l  > K2/K2  for  3 = 30°,  the_reverse  is  true  for  3 = 60°.  The 
location  of_the  maximum  value  of  K2/K2  does  not  shift,  however,  as  it 
did  for  K^/K^.  This  could  be  due  to  the  interaction  between  the  influence 
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of  the  angle  a (Fig.  1)  (which  appears  in  the  second  order  cartesian 
tensor  transformation  of  the  remote  stress  into  its  pointwise  Mode  II 
component  along  the  flaw  border)  with  the  free  boundary  effects.  (Said 
angle  a is  not  present  in  the  Mode  I transformation) . 

Effect  of  a/2m  - Tests  1 and  2 (Figs.  8,  9),  3 and  4 (Figs.  10,  7),  and 
9 and  10  (Figs.  15,  16)  represent  shallow,  moderately  deep  and  deep 
flaws  of  about  the  same  a/T!  in  pairs,  but  involving  significant  changes 
in  a/2m.  These  tests  do  not  reveal  any  identifiable  trends  insofar  as 
increasing  or  decreasing  of  the  normalized  SIF  values  with  Increasing 
a/2m.  The  most  significant  result  appears  to  be  that  K^/K^  does  not 
have  a maximum  value  at  the  point  of  maximum  flaw  penetration  in  any  of 
these  tests  but  rather  at  a/amax  of  1/3  to  1/2.  The  absence  of  trends 
in  the  normalized  SIF  magnitudes  again  suggests  a strong  role  of  the 
orientation  of  the  external  plate  surfaces  which  is  not  accounted  for  in 
the  current  normalizing  procedure. 

Effect  of  Plate  Width  - Four  tests  (13-16,  Figs.  20-22)  were  conducted 
in  order  to  assess  the  gross  influence  of  plate  width.  As  seen  in 
Figure  22,  Tests  15  and  16  indicate  the  increase  in  2m/w  from  0.31  to 
0.36  results  in  a reduction  in  K^/a  (Tra)  1/2  of  about  10%  at  the  point  of 
maximum  flaw  penetration.  However,  if  one  compares  the  average  of  the 
results  of  tests  9 and  10  (Fig.  17)  with  the  result  of  test  14  (Fig. 

21),  one  finds  the  width  influence  of  2m/w  = 0.44  to  still  be  small. 
Moreover,  comparisons  between  tests  9 and  13  (Figs.  15,  20)  (which  have 
a/Tf  and  a/2m  values  in  the  same  range  suggests  that  the  width  effect, 
while  significant  at  2m/w  = 0.64  for  K^/K^,  is  negligible  for  K2/K2. 

Photoelastic  analysis  of  secondary  slices  for  K3  determination 
revealed  that,  for  the  range  of  geometries  studied  in  the  present  pro- 
gram, the  magnitude  of  K3  was  so  small  that  it  was  swamped  by  experi- 
mental scatter  and  thus  no  K3  values  are  reported  herein. 
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SECTION  V 


DISCUSSION 

The  authors  are  aware  of  only  two  analytical  investigations  which 
appear  to  bear  in  any  way  upon  the  studies  conducted  here.  One  such 
study  (34)  conducted  by  F.  W.  Smith  and  Thresher,  focused  upon  the  part 
circular  flaw  shape  under  Mode  I loading  and  was  generally  verified  by 
the  authors  experimentally  in  an  earlier  study  (23)  for  moderate  to  deep 
flaws.  The  exact  location  of  maximum  Mode  I values  was  not  clear  from 
the  verification  tests,  however,  due  to  slight  loading  dysymmetry  in 
the  experiments.  Tests  15  and  16,  however,  (Fig.  22)  clarify  this 
point.  Using  quadratic  interpolation  of  the  results  of  (34)  they  are 
compared  with  the  Mode  I results  of  this  study  in  Figures  23  and  24. 
Differences  in  these  comparisons  would  be  attributed  to  the  difference 
in  orientation  of  the  plate  surfaces  to  the  Mode  I direction  (i.e. 
normal  to  crack  surface) . As  expected  these  differences  are  larger  for 
3 = 30°,  Fig.  23,  than  for  3 = 60°,  Fig.  24.  A second  difference  is 
that  the  analytical  study  does  not  predict  maximum  SIF  values  at  loca- 
tions other  than  at  maximum  flaw  penetration.  On  the  other  hand,  a 
second  analytical  study  by  F.  W.  Smith  and  Sorensen  (31),  noted  earlier, 
was  carried  out  on  semi-elliptical  shaped  surface  flaws  using  uniformly 
distributed  stress  distributions  normal  to  the  flaw  surface  and  parallel 
to  the  surface  dimension  of  the  flaw.  Here,  as  in  (34)  the  plate  sur- 
face was  normal  to  the  crack  surface.  These  studies  did  show  a shift  in 
the  location  of  the  maximum  normalized  SIF  away  from  the  point  of  maximum 
flaw  penetration  for  low  aspect  ratio  shallow  to  moderately  deep  flaws 
but  not  for  deep  flaws. 

In  separate  investigations  (26-30)  the  authors  have  found  that,  in 
general,  stable,  or  fatigue  flaw  growth  in  three  dimensional  problems  is 
generally  not  self-similar,  nor  do  real  flaw  shapes  assume  the  form  of 
simple  curves.  Moreover,  Sommer  et  al  (38)  have  shown  that  surface 
flaws  in  metals,  if  initiated  by  part  circular  slots,  tend  to  grow  into 
nearly  semi-circular  shapes,  then  semi-ellipses  and  finally  a bulge 
appears  near  the  region  of  (Ki/K^)  maximum  in  the  present  study  which 
distorts  the  shape  from  a semi-ellipse.  Collectively  the  above  results 
suggest  that  natural  flaw  shapes  in  three  dimensional  cracked  body 
problems  are  strongly  influenced  by  a variation  in  crack  growth  resis- 
tance along  the  flaw  border,  which  the  authors  conjecture  to  be  primarily 
due  to  a variation  in  triaxial  constraint  along  the  flaw  border.  More- 
over, the  authors  have  found  (26-30)  that,  in  complex  three  dimensional 
problems,  the  flaws  will  assume  shapes  which  tend  to  minimize  the  SIF 
gradient  along  such  borders.  However,  due  to  constraint  variations 
along  the  flaw  borders,  the  SIF  will  not  be  uniform  even  for  naturally 
grown  flaws. 

The  present  study  focused  upon  part  circular  slots  exclusively,  and 
while  such  shapes  are  commonly  employed  as  fatigue  crack  starters,  they 
cannot  represent  the  myriad  of  flaw  shapes  seen  in  real  fatigue  cracks. 
Nevertheless,  the  shift  of  the  point  of  maximum  SIF  away  from  the  point 
of  maximum  flaw  penetration  in  the  part  circular  cracks  can  be  partially 
explained  by  noting  that  the  natural  flaw  shape  is  relatively  deeper 
than  part  circular  cracks  in  the  region  where  the  SIF’s  peak  in  the  part 
circular  flaw  results.  If  such  part  circular  cracks  were  permitted  to 
grow  naturally,  they  will  tend  to  grow  more  in  such  regions  reducing  the 
local  SIF  there. 

It  should  be  emphasized  that  the  above  discussion  is  confined  to 
stable  flaw  growth  as  experienced  in  fatigue. 
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SECTION  VI 


SUMMARY  and  CONCLUSIONS 

After  developing  appropriate  analytical  expressions  for  mixed  mode 
loadings  in  three  dimensional  cracked  body  problems,  a series  of  ’’frozen 
stress"  photoelastic  experiments  were  conducted  on  plates  containing 
narrow,  part-circular  slots  inclined  at  an  angle  3 to  the  plate  surface 
and  loaded  in  uniform  simple  tension  parallel  to  the  lateral  plate 
boundaries.  Frozen  stress  data  were  extracted  from  slices  taken  mutually 
orthogonal  to  the  flaw  surfaces  and  flaw  borders,  and  inserted  into  a 
digital  computer  program  for  estimating  SIF  distributions  along  the  flaw 
borders.  These  results  were  evaluated  and  compared  with  a modified  form 
of  an  existing  Mode  I analysis.  The  following  conclusions  are  drawn: 

1)  Effect  of  3 (Fig.  1) : Although  the  influence  of  3,  or,  more  pre- 
cisely, the  angle  between  the  direction  of  the  remote  stress  and  the 
crack  plane  as  sensed  through  an  ordinary  second  order  cartesian  tensor 
transformation  at  the  crack  tip  of  the  remote  stress  was  normalized  out, 
the  effect  of  the  orientation  of  the  plate  surfaces  relative  to  the 
crack  plane  remained  in  the  data  and  significantly  altered  the  SIF 
distributions_f or  the  two  values  of  3 studied.  The  main  effect  on  the 
Kf/Ki  and  K2/K2  magnitudes  was  to  reverse  their  order  in  going  from  3 = 
30°  to  3 = 60°. 

2)  Effect  of  a/T ’ : As  a/T*  was  increased  from  0.2  to  0.7,  the  values 
of  K3/K3  were  found  to  increase  also.  Moreover^  the  SIF  distribution 
also  changed,  shifting  the  point  of  maximum  K^/Kj  further  from  the  point 
of  maximum  f law_penetration  with  increasing  depth.  Such  shifts  were  not 
observed  in  K2/K2. 

3)  Effect  of  a/2m:  Test  results  were  confined  to  0.10  < a/2m  < 0.30. 
In  this  range,  doubling  the  value  of  a/2m  did  not  produce  any  significant 
data  trends. 

4 ) Effect  of  plate  width:  Tests  were  conducted  here  and  elsewhere 
(23)  which  established  that  for  values  of  2m/w  < 0.30,  influence  on  the 
SIF  distributions  of  the  side  boundaries  of  the  plate  were  within  ex- 
perimental scatter. 

3 ) SIF  Distributions:  The  shift  in  the  point  of  maximum  K^/K^ 
away  from  the  point  of  maximum  flaw  penetration  is  a real  and  consis- 
tently observed  effect. 

6)  Role  of  K3:  For  the  geometrical  ranges  studied  here,  the  role  of 
K3  was  insignificant. 

In  applying  the  results  of  this  study,  it  is  important  to  note  that 
this  study  showed  that,  in  applying  the  frozen  stress  method  in  mixed 
mode  analysis,  where  0°  was  measured,  (rather  than  assumed)  scatter  in 
the  data  approximately  doubled  (i.e.  + 12%  rather  than  + 5%).  Moreover, 
while  these  results  are  believed  1.0  be  correct  trenas  ai  d distributions 
to  within  the  above  accuracy,  it  must  be  kept  in  mind  that  the  part 
circular  flaw  shape  is  not  necessarily  the  preferred  shape  of  a naturally 
growing  fatigue  crack,  and  consequently,  may  lead  to  higher  SIF  gradients 
along  the  flaw  borders  than  for  other  flaw  shapes.  Nevertheless,  these 
results  are  believed  to  be  the  most  practical  and  comprehensive  for  the 
problem  studied  that  are  currently  available.  The  authors,  however,  do 
not  recommend  extrapolation  to  significantly  different  flaw  geometries. 
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Figure  1)  Problem  Geometry 
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Figure  4)  Determination  of  K*  from  Typical  Slice  Data 
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Figure  5)  Kelvin  Model 
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Figure  6)  Typcial  a)  Mixed  and  b)  Mode  I fringe  patterns 
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Table  I 


T' 


T 

sin  8 


2 1/2 
* a sin  8(Tra; 

1/2 

* a sin  8 cos  8(ira) 


Test 

No. 

6 

a (mm) 

2R(mm) 

2m  (mm) 

T(mm) 

T* (mm) 

a/2m 

a/T’ 

a e 
max 

Notch.  . 

o(kpa) 

1 

30° 

5.49 

76.2 

39.37 

13.67 

27.43 

0.139 

0.200 

31.11 

0.889 

103.10 

2 

30° 

6.09 

38.1 

27.94 

13.72 

27.43 

0.218 

0.222 

47.17 

0.533 

102.30 

3 

60° 

4.32 

101.6 

40.99 

13.46 

15.54 

0.105 

0.278 

23.79 

1.372 

76.76 

4* 

60° 

4.32 

25.4 

19.08 

12.78 

14.76 

0.226 

0.293 

48.70 

0.432 

80.28 

5 

30° 

9.53 

101.6 

59.18 

13.59 

27.18 

0.161 

0.350 

35.66 

1.372 

93.72 

6 

30° 

10.89 

76.2 

53.34 

13.46 

26.90 

0.204 

0.405 

44.43 

0.889 

104.10 

7 

60° 

7.26 

50.8 

35.56 

12.75 

14.73 

0.204 

0.493 

44.43 

0.940 

81.03 

8 

60° 

7.62 

101.6 

53.52 

13.64 

15.75 

0.142 

0.484 

31.79 

1.372 

75.72 

9 

30° 

9.53 

76.2 

50.29 

7.32 

14.63 

0.189 

0.650 

41.40 

0.889 

63.45 

10 

30° 

9.53 

50.8 

33.05 

7.29 

14.58 

0.288 

0.653 

51.31 

0.940 

66.07 

11 

60° 

10.99 

76.2 

53.57 

14.17 

16.38 

0.205 

0.672 

44.66 

0.889 

72.28 

12 

60° 

10.16 

50.8 

40.64 

12.45 

14.38 

0.250 

0.707 

53.13 

0.940 

83.03 

13 

30° 

17.60 

203.2 

114.30 

14.05 

28.09 

0.154 

0.627 

34.23 

3.429 

72.90 

14 

30° 

18.69 

101.6 

78.74 

14.30 

28.60 

0.237 

0.654 

50.81 

1.372 

71.59 

15a 

90° 

6.35 

76.2 

42.11 

12.70 

12.70 

0.151 

0.500 

33.56 

0.889 

66.10 

16b 

90° 

6.35 

76.2 

42.11 

12.45 

12.45 

0.151 

0.510 

33.56 

0.889 

67.74 

* Avg.  of  4 experiments 


All  plates  are  178  mm  wide  except 
a - plate  width  - 137.8  mm 
b - plate  width  ■ 117.2  mm 


Test  Geometries  for  Figures  7-22 
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a/2m  = 0.226 


J 1 1 i i i 

1/6  2/6  3/6  4/6  5/6  I 


ax 


Figure  7)  Mixed  Mode  SIF  Distributions  for  Tests  4 
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p-Ip* 


Figure  8)  Mixed  Mode  SIF  Distributions  for  Test  1 
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P-lP- 


1.0 


Test  2 
see  Table  I 

$ = 30° 
a/T'=  0.222 
a/2ms0.2l8 


Figure  9)  Mixed  Mode  SIF  Distributions  for  Test  2 
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P-lP- 


^max 

Figure  10)  Mixed  Mode  Distributions  for  Test  3 
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Figure  11) 


Mixed  Mode  STF  Distributions  for  Test 


mu  a 
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P-|P- 


1.0 


0 = 30° 
a/T'=  0.405 
a/2m  = 0.204 


i i i i l I 

1/6  2/6  3/6  4/6  5/6  I 


Figure  12)  Mi::ed  Mode  SIF  Distributions  for  Test  6 
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P-lP- 


Figure  13) 


Mixed  Mode  SIF  Distributions  for  Test 


I I IUA 
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Figure  14)  Mixed  Mode  SIF  Distributions  for  Test  8 
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o I 1 1 I 1 l 

1/6  2/6  3/6  4/6  5/6  I 


Figure  15)  Mixed  Mode  SIF  Distributions  for  Test  9 
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Figure  16)  Mixed  Mode  SIF  Distributions  for  Test  10 
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P-lP- 


Figure  17)  Mixed  Mode  SIF  Distributions  for  Avg. ^ Tests  9,  10 
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1.0 


0.5 


K2/K2 


Test  1 1 
see  Table  I 
0 = 60° 


a/T'  = 0.672 
a/2  m = 0.205 


i i i l I — 

1/6  2/6  3/6  4/6  5/6 


Figure  18)  Mixed  Mode  SIF  Distributions  for  Test  11 
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P-|P- 


Figure  19)  Mixed  Mode  SIF  Distributions  for  Test  12 
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p-Ip- 


Figure  20)  Mixed  Mode  SIF  Distributions  for  Test  13 
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Figure  21)  Mixed  Mode  SIF  Distributions  for  Test  14 
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F~  |P" 


K,  /cr(77*a) 


Figure  22)  Mixed  Mode  SIF  Distributions  for  Test  15,  16 
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P-F 


K,/K, 


MTS 


Figure  23)  Comparison  of  Mode  I Results  with  Modified  Thresher-Smith  (MTS) 
Theory  (6  = 30°) 
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K|/K, 


Figure  24)  Comparison  of  Mode  I Results  with  Modified  Thresher-Smith  (MTS) 
Theory  (g  = 60°) 
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